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Abstract. We present an algorithm for efficient computation of the constant 
term of a power of a multivariate Laurent polynomial. The algorithm is based on 
univariate interpolation, does not require the storage of intermediate data and 
can be easily parallelized. As an application we compute the power series expan- 
sion of the principal period of some toric Calabi-Yau varieties and find previously 
unknown differential operators of Calabi-Yau type. 



1. Introduction 

The constant terms of powers of Laurent polynomials appear in different areas 
of mathematics and theoretical physics I'D | , | DvK| . They show up prominently in 
the expansion of the fundamental period 



If 1 dXi dX 2 dX N 



(2iri) N J r 1 - zf(X) X 1 X 2 X 



N 



of a Laurent polynomial in X\, . . . , X N , and r is the TV-dimensional torus cycle 
\Xi\ — e. In the theory of mirror symmetry for Calabi-Yau hypersurfaces in a toric 
variety, this fundamental period is of great importance, IBKI , IBvSI . To compute 
the Taylor expansion of <f> in the neighborhood of z = we write 

1 ^(/W)v 



1 - zf(X) ^ 

and integrate using Cauchy's formula. The n-th coefficient in the power series 
expansion of <i> is therefore the constant term of the n-th power of the Laurent 
polynomial /. For example the constant term of the 151st power of 

X ZT 

f = 1 h ZT + T -\ 

J YZ XY 

1 X ZT 

+ Z +t+X+ y + — 

is 

15412036066982883611159466717890839926274227993361685769096965357956125 
0836097113850549748895583119569242295079022614473032754474202469738117581 
03097074502829198076370950235391810731785760778732696320, a number with 200 
digits. 

In this paper we consider the following general problem: 
How to compute for a multivariate (Laurent) polynomial f the coefficient of the 
monomial X a , a = («i, . . . , a n ) ofthep-th power, p e N off. 
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At first glance this problem is trivial. Why cant't we just split the problem in 
smaller ones and compute f r =: g = J2 bpXP and f s ='-h = Y^ c^X^ s.t. r + s = p 
and get the coefficient we are interested in by Cauchys product formula 

(1) coeff. ofX a inf p = ^ b iCj . 

A general polynomial in n variables in which the degree of each individual vari- 
able is not greater than d has (d + 1)" monomials. In our example n = 4, p = 
151, d = 2 after multiplication with the common denominator XYZT. In or- 
der to use (I| we could compute f 75 and f 76 which will require to store about 
2(2 • 76 + 10 9 coefficients of intermediate polynomials g and h. While this 
number is not really big, we must keep in mind that, in general, the absolute 
value of the coefficients grow exponentially as p increases. Therefore, the amount 
of memory needed to store one particular coefficient will, typically, linearly grow 
with p. Another possibility is to use a multivariate discrete Fourier transform: 
evaluate the polynomial / on sufficently many points of the n-dimensional grid 
of roots of unity of sufficient high degree, compute the p-th power of the values 
and apply the inverse discrete Fourier transform to get the single coefficient we 
are interested in. Here, again, the number of interpolation nodes and therefore 
of the intermediate values which have to be stored is of the same order as the 
number of monomials in f p i.e. (2 • 151) 4 . 

In this paper, we present an efficient and, we believe, simple algorithm which 
adresses this problem. Our approach is similar to the DFT method in that it 
based on evaluation and iterpolation. However, the amount of memory needed 
is much lower and the running speed can be improved by a trivial parallelization 
of the evaluation. We also describe the implementation of the parallel version of 
the algorithm for the CUDA-enabled graphics hardware. Finally, an application 
to the problem of finding the Picard-Fuchs differential equation of the family of 
toric Calabi-Yau varieties is given. This problem was the main motivation for this 
work. 



Acknowledgement: I like to thank D. van Straten for suggesting the problem 
and E. Schomer for useful discussions. This work was realised with support from 
the DFG priority program Algorithmic and Experimental Methods in Algebra, Ge- 
ometry and Number Theory, as part of the project Monodromy algorithms in Sin- 
gular. 

2. Algorithm 

2.1. Notations and Preliminaries. Let / be a multivariate polynomial in vari- 
ables Xi,...,Xk with rational coefficients. Using multiindex notation we write 

iGA 

i := (ii, . . . , i k ) and AcZ'a finite index set. Let deg,,(/) denote the degree of / 
in the variable X r . For a given fc-dimensional multiindex j we denote by [/] . the 
coefficient of / in front of X 3 i.e. 
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Furthermore let f denote the (fc— 1) -dimensional multindex obtained by trun- 
cating j after the (fc — l)-th component and for an r-dimensional multiindex i let 
ij denote the concatenation of both i.e. for j = (ji, . . . ,ju) audi = (ii, . . . , i r ) 

f ■ = O'i, • ■ ■ , jk-i), 

ij ■ = (ii,...,i r ,ji,...,jk)- 

Given N + 1 pairs (u a , w a ) e Q 2 , a = 0, . . . , N such that u a ^ up for a ^ (i, 
there is a unique univariate interpolation polynomial of degree N such that 

Ijv(« Q ) = w a , Va G {0, . . . ,N} 

N 

1 n (uq, u N ; io , ... , w N )(Y) = ^2 Wih(Y), 



i=0 



where 



is the Lagrange basis polynomial for v*. In has rational coeffiients andl N (u a ) 
w a , Va = 0, . . . , N. 

Let f u denote the polynomial / with u substituted for the last variable i.e. 



fu{Xi, . . . , X n —x) — f{X\, . . . , X n _i, u). 

We have the following observation about the coefficients of the interpolation 
polynomial and those of /„. 

Lemma 1. Let f e Q[X U . . . , X n ], N := deg n (f), {(u a , w a )} a=Q ^ N as above, and 
i = (ii, . . . , i n ) e N™. Then, we have 

(2) [f]i= [l N {uo,---,u N ;[f ao ]i',...,[f UN ]i')] in . 
Proof. Let 

(3) f(Xi,...,X n )= ^ Cviv*...v n X'i~ ■ ■ • X„ n , 

U\,V^,...,V n 

and 

JV 

(4) l N (u ,...,u N ;[f Uo ] i ,,...,[f UN ] il )(Y)=J2akY k . 

k=0 

Then, the right hand side of {2j equals a ln and, therefore, we have to show that 
= a in . The coefficients ao, . . . , a n of the interpolation polynomial are found 
by multiplying the vector ([/„„],/, [f Ul }i',- ■ ■ , [fu N ]i>) T by the inverse of the Van- 
dermonde matrix for ti ,...,%. This inverse exists since u , . . . , un are pairwise 
distinct. Therefore, we have 



(5) 





ao 
a i 


\ 




( 


1 
1 


Uq 
Ul 


"0 

uf 


Uq 
■■■ < 


\ 


-1 


( 


[fu ]i' 
[/"ill' 


\ 


V 


a N 


J 






1 


UN 




• • • U N 

N 


J 






[fu N ]i' 


J 



4 



PAVEL METELITSYN 



By l|3} we have 



[fuj ]i' 



T U j JL 1 A 2 . 



vi,V2,...,u n —i,r 
Ciii 2 ...i,i-iTMj 

T 

^ ] c i'T u i , 



where in the last expression i'r is an abbriviration for the multiindex (ii , 
Using this, rewrite {5} as 



(6) 
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At the same time 
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and with {3} 
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follows. 

Since deg n (/ p ) = p ■ deg„(/) we have 

(8) m = [l w («o, • ■ • ,u N ; [/£>, ... , [fUi')] in > 

where AT = p-deg„(/). 

2.2. The algorithm and its application to the constant term series of a Lau- 
rent polynomial. Let h e Q[Xf , . . . , X^] be a Laurent polynomial. We want 
to calculate the first coefficients of its constant terms series ([/i p ] ) p eN- Let S = 

f Sl X Sk s- 

h. Then S7i G 



X k ] and 

where p ■ s := (psi, . . . ,ps„). We want to apply {8} to / = Sh and i = p ■ s. 

The simple idea now is to compute f£. , j = 0, . . . , N which appear in |8 
recursively i.e. 



(9) 



J U i 



1 (uo,...,u N ;[f* 
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and so on for f UjUz , f Uj u z u k , With each level of recursion the number of vari- 
ables decreases by one, therefore the depth of the recursion equals n. For sim- 
plicity assume 

de gl (/)=---=deg„ (/) = <! 
The computation of [f p ] i proceeds as follows: 

(1) For N = dp fix u , ■ • • , u N as in Lemma and compute V^ 1 = (T4,-)> where 
V is the Vandermonde matrix of (u , . . . , u N ). 

(2) Invoke the recursive procedure COEFF(z, p, k, A) with initial parameters 
k — n and A ^"coefficient matrix of /". 

The pseudo code of the procedure COEFF is given in the Algorithm 1. 

Algorithm 1 COEFF(fc, i, A,p) 

Input: k > 0, i e N k , A e Matj d x ■■■ x d) , p>l 

Output: [/p], 
i: for all < s < N do 

2: />', <- Do it using Horner's rule 

3: if k = 1 then 

4: w s <- B' p 1 1 B has k — 1 = indices i.e. B is a scalar 
5: else 

6: w s <- COEFF(fc — 1, B,p) 
7: end if 
8: end for 

9: S <- J2f=o w sV~i k s II Compute the «fc-th component of V~ 1 w 
10: return S 



To improve the performane we apply some well known computational tricks. 
The evaluation in the line 2 is done using the Horner's rule which for a polyno- 
mial h(x) = ao + a\x + • • ■ + aix 1 at x = xq says 

f(x ) = a + x (ai + x (a 2 + x (. ■ ■ )))■ 

In this way the expensive computation of powers of x is avoided. The Vander- 
monde matrix V can be inverted using classical Gauss elimination algorithm in 
0(N 3 ) time. Since on each level of recursion we need only one row of V^ 1 (line 9 
in the code above) precomputing the whole inverse seems too expensive. Fortu- 
nately if the nodes uq,...,un are choosen to be equidistant there is a recursive 
algorithm which computes one single row of V^ 1 mO(N 2 ) (T). The computation 
of the p-th power in the line 4 can be done using binary powering algorithm as 
described in (K)- We also note that w s 's in line 6 do not depend on each other 
and thus the order in which they are computed does not matter. In particular 
they can be computed parallely. Finally the number of nodes TV need not to be 
the same on each recursion level but depend on the degree of / in a particular 
variable i.e. we have to choose N k ■= max{deg 1 (/), . . . ,deg fc (/)}-ponthe (n—k)- 
th recursion level. For example if / has degree 2 in one variable and degree 1 in 
the remaining variables we could improve speed by considering 2p interpolation 
nodes only on one level of recursion and p interpolation nodes on the ramaining. 
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If / happens to have special form 
(10) f = A+ BXx + CX\ 

with A, B, C e Z[X 2 , • ■ • , X k ] a trick can be applied to slightly speed up the com- 
putation. If 1 10 1 holds, then so is for /„, 

fu = A u + B u Xi + C U X\ 

with A u , B U ,C U € Z[X 2 , . . . , Xk-i]. Consequently, this decomposition holds on 
all levels of recursion. We note that in this case 



and 
(11) 
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X{ i.e. 



j. We apply this "short- 
when we have 



since we need only the terms with X\ t+j 
cut" on the last but one level of recursion to compute v 

to deal with polynomials in two variables. We use the modified version of COEFF 
and a new procedure SPLIT2 for this. Thus the depth of recursion is reduced by 



one. For the case n = 4 the call tree of the algorithm is depicted in Figure 2.2 



Algorithm 3 Modified version of COEFF 
Input: k > 0, i € N k , A € Matj dx ■ ■ ■ X d) , p>l 

Output: [f% 
i: for all < s < N do 

2 : Bii,.. -,ih-i ^~ Ai l ^... t i k _ l .jU 1 s 

3: if k = 2 then 

4: w s <r- SPLIT2(i',B,p) 
5: else 

6: w s COEFF(fc — 1, i',B,p) 
7: end if 
8: end for 
9: S <- HfLo w sV is 
10: return S 



The bottom procedure SPLIT2 is called exactely TV" 2 times. Inside the SPLIT2 
procedure we have two nested loops of depth 2. The bodies of both inner loops 
(lines 12-14 and 25-28 in the listing ij are executed both at most ^ times. There- 
fore the time complexity of the algorithm is of order 0{N n ). Thanks to SPLIT2 we 
gain a factor ^^=2 compared to the version which does not use SPLIT2. If we take 
into consideration the linear growth of the length of the numbers involved into 
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Algorithm 4 SPLIT2(i, A,p) 

Input: k > 0, i e N fc , A e Mat(d x •• • x d) ,p> 1 

k 

Output: [fP]. 

1: s <- 

2: m <- f + 1 

3: for all < i < AT do 

4: <- (A 00 + Ui(A 01 + UiA 02 )) ■ (A 20 + u l (A 21 + UiA 22 )) 

5: t <- A w +Ui(Au +UiA 12 ) 

6: if n = mod 2 then 

7: Um-l,i <- 1 

8: else 

9: U m -i,j <- i 

10: end if 

11: t <- t 2 

12: for all < j < to - 2 do 

13: V rn — 2 —j^i i t • Vrn—l—j.i 

14: end for 

15: w'i 4- 1 

16: w 0iJ 

17: end for 

18: t <r- 

19: for all < i < N do 

20: S i— S + CiZi 

21: end for 

22: s <- M ■ s 

23: for all 1 < j < m do 

24: t <- 

25: for all < i < N do 

26: w't 4- w't ■ Wi 

27: t 4— t + ■ Wi ■ Ci ■ Vj t i 

28: end for 

29: S <- S + Mj ■ t 

30: end for 

31: return s 



computation we arrive at assymtotical running timepiV™ or d n p n+1 , since N de- 
pends linearly on p and d. 

As we can see from the description of the algorithm, the memory consump- 
tion is O(N) = 0(dp), which is the space needed to store the n - 1 rows of the 
inverse Vandermonde matrix. In particular, for fixed d it is independent from the 
number of monomials / has. 

Remark 1. When d is fixed, the running time depends on p and the number of 
monomials in /. 

2.3. Implemetation details. As was mentioned before, the algorithm can be par- 
allelised in a very straight-forward way. As a target platform we choose a graphics 
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Figure 1 . Call tree of the algorithm for k = 4. 



card built around NVIDIA GeForce GTX 480 GPlf] This GPU is capable of run- 
ning up to 480 parallel execution threads. The GPU is programmed using an 
extension of C/C++ programming language called NVIDA CUDA^ C |N). A typi- 
cal CUDA program is executed on both the CPU and the GPU. The parts of the 
code that run on the graphics hardware are called kernels and each kernel can be 
run TV times in parallel by N different CUDA threads. The threads are organized 
in one-, two- or three-dimensional thread blocks while blocks are organized in 
one- or two-dimensional grid. The graphics hardware was origianlly created for 
the purpose of low-precision floating point operations it uveils its full comput- 
ing power only when dealing with 32-bit floating point numbers. Therefore we 
can perform exact computations only on numbers within the range — 2 23 . . . 2 23 
or equivalently -8388608 . . . 8388608. That is because only 23 bits of 32 are used 
to store the significant digits (8-bits being reserved for the exponent and one re- 
maining bit for the sign). This range is far too small to be useful due to the rapid 
growth of the numbers [f p ] . Therefore we used a residue number system (RNS) to 
represent large integers. Given a (fixed) set of pairwise coprime natural numbers 
mi,...,m s , the integer x < M := J] m l is represented by the system of residues 



xi : 



x mod mi 



mod 



M is called the dynamic range of the system. By the well-known Chinese Re- 
minder Theorem the number x can be reconstructed from (x%, . . . ,x s ; mi, . . . , m s ). 
In practice, the conversion from an RNS representation of a; to a decimal repre- 
sentation is done with the Mixed Radix Conversion (MRC) i.e. x is represented in 
the form 



di + a 2 TOi + a 3 mim 2 



a s m\rri2 ■ 



Graphics Processing Unit 
2 Compute Unified Device Architecture 
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Figure 2. Arrangement of threads. 



where en's are called mixed radix digits. Once a/s are known the decimal value 
can be computed easily. One classical algorithm for finding en's from RNS repre- 
sentation of x is given in |SR) . 

In our implementation the program branches in several execution threads in 
the procedure COEFF on the top level of recursion. The threads are organized in 
b blocks with t threads per block. All the threads within the fc-th block performs 
computation modulo prime number p k . The i-th thread of fc-th block computes 
w-j mod p k for j = i mod t. The arrangement of threads is illustrated in Figure 
2. 



3. Application to toric Calabi-Yau varieties 

In the theory of toric Calabi-Yau varieties the constant terms of powers of Lau- 
rent polynomials occur in the following way. Consider a reflexive lattice polytope 
A c M. 4 . Every vertex Uj = (v^i, . . . , v it4 ) e Z 4 of A corresponds to a monomial 
n*=i X] u ■ Then the Laurent polynomial corresponding to A is 

i=i 3=1 

This polynomial defines a family of hypersurfaces := {1 - z ■ /a = 0} c (C x ) 4 . 
The function 

oo 

(12) 

»=0 

is the period of the holomorphic 3-form 

( 1 dX 1 A dX 2 dX 3 dXi 



l-zf A (X) X 1 X 2 X 3 X 4 
on X z and is a solution of the Picard-Fuchs equation 

L$(z) = 
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where 

(13) L = P o (0) + zPx{6) + ■ ■ ■ + z k P k {9), 

6 := z-j; and Pi are polynomials with integral coefficients 

di 
3=0 



The numbers (a,),^ 
(14) 



'[/*]„> *>0 

0, otherwise 



then satisfy the following recurrence relation 

(15) Po(n)a n + Pi{n - l)a„_i + • • • + P fc (n - k)a n - k = 0, 

for all n e N . We say that the recurrence has length k + 1 and degree d := 
max {deg Pi}. Conversely, if a sequence satisfies (15) then its generating function 
is anihilated by the differential operator (13). 

Given A we want to find the corresponding differential operator L. To do this, 
we compute first N coefficients of the constant terms series (12) of /a and try 
to find a reccurence of the form (15). Assume there is one with lenght k + 1 and 
degree d. Then it is determined by (k + 1) (d + 1) coefficients b LJ , i = 0, . . . , k, j = 
0, . . . , d of Pq, . . . , Pk- We check if there is a solution to system of TV + 1 linear 
equations for the (k + l)(d + 1) unknowns b l3 

P (N)a N + P 1 (N-l)a N _ 1 + --- + P k {N-k)a N - k = 



(16) 



P (2)a 2 + Pi(l)ai + P 2 (0)a = 

P o (l)ai+Pi(0)a o = 

P o (0)a o = 



We have to keep this system overdetermined i.e. to have 7V + 1 > (fc + l)(d + l). If 
there is a solution which does not change when we add more equations, we hope 
that we found the correct operator. Since, a priori, we do not know the length 
and the degree of the reccurence we make the Ansatz (16) for all k, n which are 
allowed by the requirement that (16) be overdetermined. Thus it is essential to 
be able to compute [/'] for as many i's as possible in reasonable time. 

This method was used by Batyrev and van Straten IBvSI and has been polu- 
lar ever since. More recently, Batyrev and Kreuzer used this approach in IBKI 
to determine the Picard-Fuchs operators for some new families of Calabi-Yau 
threefolds with Picard number 1. They succeeded in 28 cases out of 68. Using 
our algorithm were able to compute 4 more Calabi-Yau operators. Although we 
knew up to 300 constant terms, in many cases the recurrence did not show up. 
The operators and their corresponding polynomials are: 

#24: 

f — I^V-lZ-l^Tj- zt_ , i . x , y , x , xy_ _r_ y_ i t_ _i_ _t_ i y_ i 1_ i 

724 T' J T^-|-jf-t- XY T z T z T ZT T ZT T ZT T T T y T X Y ' X ' X ' 
+ Y ^ ~T ~ZT Z -\- y 
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D 24 = 97 2 8 4 + 97z9 (-291 - 1300 9 - 2018 9 2 + 1727 9 3 ) 

+z 2 (-2709792 - 10216234 6*- 16174393 6> 2 - 13428812 9 3 - 1652135 8 4 ) 
+z 3 (-138000348 - 443115594 8 - 568639497 9 2 - 364126194 6> 3 - 81753435 6 / 
+z 4 (-3049275024 - 8869415520 9 - 10006378570 9 2 - 5423394464 9 s - 1175502862 9 A ) 
+z 5 (-38537290992 - 103964102350 9 - 106108023451 9 2 - 50507429234 9 3 - 9726250397 9 A ) 
+z 6 (-308040167808 - 781527778884 9 - 733053660150 9 2 - 312374434824 9 3 - 52762935894 9 4 ) 
+z 7 (-1619360309088 - 3901093356168 9 - 3399527062044 9 2 - 1313199235080 8 3 - 195453433908 6> 4 ) 
- 144z 8 (6» + l) (3432647479 8 3 + 22487363787 9 2 + 5080861471 1 9 + 38959393614) 
-432z 9 (9 + 2) (9 + 1) (1903493629 9 2 + 10262864555 9 + 14314039440) 
-438048z 10 (186298761 + 5992902) (9 + 3) (8 + 2) (9 + 1) 
-368028363456Z 11 (9 + 1) (8 + 2) (9 + 3) (6 + 4) 

#39: 

/39 = yz + §-y- + ZT + T + jp + y + jff + ^ + + y+ t+ Z+ t+ X+ y + ^ + yt 

~ ZT ~ X ~ Y ~ X ' Y ~ Z 

D 39 = 16 9 4 -4z9 (12 + 53 9 + 829 2 + 29 3 ) 

+z 2 (-5120 - 18308 9 - 26199 9 2 - 18410 9 3 - 4895 6> 4 ) 

+z 3 (-143808 - 430092 6 1 - 497452 2 - 272424 9 3 - 60679 6» 4 ) 

+« 4 (-1478544 - 3987101 9 - 4034628<9 2 - 18708 3 8 6I 3 - 344527 6» 4 ) 

-z 5 (9 + 1) (1076509 9 3 + 5847783 9 2 + 11226106 8 + 7492832) 

-2z 6 (8 + 2) (9 + 1) (944887 9 2 + 4249317 9 + 5045304) 

-3328z 7 (518 9 + 1381) (9 + 3) (8 + 2) (9 + 1) 

-621920z 8 (8 + 1) (6» + 2) (8 + 3) (8 + 4) 



#41: 



fa — XT + YZ + ZT + + T + + + -^y + -^j, + XZT + yy + Y ZT + XY z xyt 
+ JTZT + yzt + XYZT + XYZT + Y + X + Z+ T + z + Y + ¥ 

D 4 i = 91 2 4 + 91261 (-273 - 1210 8 - 1874 9 2 + 782 9 3 ) 

+z 2 (-2649920 - 9962953 9 - 15227939 6» 2 - 11622522 9 3 - 2515785 9 4 ) 

+z 3 (-110445426 - 348819198(9 - 4326078686» 2 - 2586781266 13 - 59827597 8 4 ) 

+z 4 (-1915723890 - 5439732380 8 - 5901995820 8 2 - 2998881218 8 3 - 612043042 6» 4 ) 

+z 5 (-18479595006 - 48522700563 8 - 47503242813 6» 2 - 21226829058 9 3 - 3762840342 6» 4 ) 

+z 6 (-110147546634- 271941545379 9- 244753624741 9 2 - 98210309094 8 3 - 15265487382 8 4 ) 

+z 7 (-422269162452 - 991829482602 8 - 831965057114 6> 2 - 304487632282 6» 3 - 42103272002 6> 4 ) 

-2z 8 (8 + 1) (39253400626 8 3 + 275108963001 8 2 + 654332416678 8 + 521254338620) 

-z 9 (8 + 2) (8 + 1) (94987355417 6> 2 + 545340710193 8 + 799002779040) 

-15402 10 (43765159 8 + 149264765) (8 + 3) (8 + 2) (8 + 1) 

-2 2 35 2 7 2 ll 2 1197l2 n (6» + 1) (8 + 2) (8 + 3) (8 + 4) 



f _ X 1 ZT 171 V 4-T 4- — 4-Y 4- — 4- — 1 1 1 Y , _, 1 , 1_ , 1 , T 

J 38 y'Y' ' ' 'Z' 'X'X' XZT ' XZT ' XT ' ZT ' Z ' YZ ' Y 



#38 

/ 

D 3 s = 102 2 6» 4 - 102z (204 + 911 8 + 1414 6» 2 + 116 8°) 

+z 2 (-2663424 - 9947652 9 - 14508941 6» 2 - 9892670 9 3 - 2596259 6> 4 ) 
+2 3 (-67967496 - 206933112 9 - 239004708 9 2 - 125234088 9 3 - 25685301 6» 4 ) 
+2 4 (-598491604 - 1608054100 9 - 1587508748 9 2 - 687051032 8 3 - 112357900 8 4 ) 
+2 5 (-2495389956 - 6085656898 8 - 5273754198 6> 2 - 1927713868 8 3 - 254678692 8 4 
+z 6 (-5385015134 - 11995897911 8 - 9101625228 8 2 - 2758627602 6> 3 - 283337071 8 
+z 7 (-5612134720 - 11209872916 8 - 7075746650 8 2 - 1555791344 8 3 - 86504770 8 4 
+122 8 (6» + 1) (7613560 6» 3 + 27844427 8 2 - 51849552 6* - 134696600) 
+2 9 (8 + 2) (8 + 1) (60585089 9 2 + 495871401 8 + 595115780) 
-6OO2 10 (10279 - 113205) (8 + 3) (8 + 2) (9 + 1) 
-67900002 11 (6» + 1) (8 + 2) (9 + 3) (8 + 4) 



+ — + — + — + - + - + - + - 

~YZ ~ ZT ~ XT 'X'Y'X'T 
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